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ABSTRACT 

Methyltransferases that use S-adenosylmethionine 
(AdoMet) as a cofactor to catalyse 5-methyl uridine 
(m^U) formation in tRNAs and rPNAs are widespread 
in Bacteria and Eukaryota, and are also found in 
certain Archaea. These enzymes belong to the 
COG2265 cluster, and the Gram-negative bacterium 
Escherichia coli possesses three paralogues. These 
comprise the methyltransferases TrmA that targets 
U54 in tRNAs, RImC that modifies U747 in 23S rRNA 
and RImD that is specific for U1939 in 23S rRNA. The 
tRNAs and rRNAs of the Gram-positive bacterium 
Bacillus subtilis have the same three m^U modifica- 
tions. However, as previously shown, the m^U54 
modification in B. subtilis tRNAs is catalysed in a 
fundamentally different manner by the folate- 
dependent enzyme TrmFO, which is unrelated to 
the E. coli TrmA. Here, we show that methylation of 
U747 and U1939 in B. subtilis rRNA is catalysed by a 
single enzyme, YefA that is a COG2265 member. A 
recombinant version of YefA functions in an E. coli 
m^U-null mutant adding the same two rRNA methy- 
lations. The findings suggest that during evolution, 
GOG2265 enzymes have undergone a series of 
changes in target specificity and that YefA is closer 
to an archetypical m^U methyltransferase. To reflect 
its dual specificity, YefA is renamed RlmCD. 

INTRODUCTION 

Ribothymidine (m^U) is present in tRNAs from the three 
domains of life (1-3) and is additionally found in the 



rRNAs of all Bacteria (4; http://www.ecosal.org) and 
certain Archaea (5). The RNAs of the best-studied bacter- 
ium, Escherichia coli, contain three m^U's. Two of these 
are in 23S rRNA and one in tRNA, and these modifica- 
tions are conserved to varying degrees in other bacteria. 
The most highly conserved m^U is at position 54 in the 
T-loop of tRNA (6) and in E. coli this modification is 
added by the methyltransferase TrmA, which was 
formerly termed RumT (7). In 23S rRNA, m^U at nucleo- 
tide 1939 is present in all bacterial ribosomes characterized 
to date (8), and its formation is catalysed by the E. coli 
enzyme RlmD (formerly RumA) (9,10). The third ni^U, at 
23S rRNA position 747, is the least conserved of the three 
and is found mainly in Gram-negative Beta-, Epsilon- and 
Gamma Proteobacteria (3) and, in E. coli, this modifica- 
tion is added by the methyltransferase RlmC (formerly 
RumB) (10). All three E. coli m^U methyltransferases 
use 5-adenosyl-L-methionine (AdoMet) as the methyl 
group donor, and their sequence similarity indicates that 
they share a common evolutionary origin (3). The se- 
quences cluster within the family of orthologous groups 
COG2265 (11), which contains homologues from phylo- 
genetically distant organisms such as the tRNA m^U54 
methyltransferases Trm2p, found in the yeast 
Saccharomyces cerevisiae (12), and PAB0719 in the 
hyperthermophilic archaeon Pyrococcus abyssi (3). 

The sites of m^U modification in the RNAs of Bacillus 
subtilis are similar to those in E. coli — the tRNAs of B. 
subtilis contain m^U54, and our earlier provisional studies 
on B. subtilis 23S rRNA indicated that nucleotides U747 
and U1939 {E. coli numbering system) were also 
methylated. However, this similarity in m^U location 
obscured the fact that there are fundamental differences 
in the modification mechanisms in B. subtilis and E. coli. It 
had already been shown that in B. subtilis, the tRNA 
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m^U54 modification is added by a folate-dependent 
metliyltransferase (TrmFO) that bears no relationship to 
the COG2265 family members (13), and this discov- 
ery begged the question as to how m^U's might be 
added to the rRNA. Our searches of the B. subtilis 
genome revealed no paralogues of TrmFO, although two 
COG2265-type genes, yef'A and yfjO, were evident. The 
protein products YefA and YfjO showed sufficient 
sequence similarity to each other, and to the E. coli 
RlmC and RlmD enzymes, to rank them as the most 
likely candidates for catalysis of m^U modification in B. 
subtilis rRNA. 

In this study, we have used Matrix Assisted Laser 
Desorption/Ionization (MALDI) mass spectrometry 
(MS) to demonstrate conclusively that m^U's are present 
at nucleotides 747 and 1939 in B. .subtilis 23S rRNA. The 
rRNAs were then analysed after inactivation of yefA or 
yfjO to establish whether the protein products of these 
genes were responsible for the m^U modifications. Enzyme 
function was re-established by expressing plasmid-encoded 
recombinant protein in the B. subtilis null-mutant and also 
in an E. coli mutant in which all three trmA, rlmC and 
rlmD genes had been deleted. Our findings show that 
compared to the E. coli RlmC and RlmD counterparts, 
a functional shift has occurred during the evolution of the 
B. subtilis YefA and YfjO methyltransferases. 

MATERIALS AND METHODS 

Bacillus subtilis and E. coli strains 

Bacillus .subtilis 168 (trpC2) was used here as the reference 
strain with a wild-type rRNA modification pattern. The 
B. subtilis knockout strains, BG12826 (>'e/L4::pMutin2) 
and BG12911 (j'//0::pMutin2), originated from the collec- 
tion of the Japanese/European B. .subtilis Functional 
Analysis Program, and were created by a standard single 
crossover-based protocol using PCR-aniplified fragments 
of target genes cloned in the pMUTIN2 vector (14). The 
AyefA::Cm deletion mutant was derived from the 
B. .subtilis 168 strain. The upstream 891 bp and down- 
stream 919 bp fragments of the yefA gene were PGR 
amplified from the genomic DNA using the primers 
yef A-up 1 (5'-GCGCTCGGAATCCC ACGCG AGG) 
together with yefA-up2 (5'- CGACCTGCAGGCATGCA 
AGCT CTATAATCACCTCTAGTTTTAAGCGAA 
AG), as well as yefA-downl (5'- GCTCGAATTCACTGG 
CCGTCG TAAATAGTCAAAAATCTATGCTGATG 
G) together with yefA-down2 (5'-TAGTCCCTCGCAGC 
TTGCCCGC). The sequences underlined in the yefA-up2 
and yefA-downl primers are homologous to the ends of 
the pUC19-cloned 727-bp chloramphenicol resistance cas- 
sette (15) that was amplified using the PGR primers Cm5 
(5'-CGACGGCCAGTGAATTCGAGC) together with 
Cm3 (5'-AGCTTGCATGCCTGCAGGTCG). The three 
PGR products were mixed in equimolecular amounts and 
joined in a subsequent PGR. The resulting 2.5-kb DNA 
fragment was used to transform competent B. .subtilis 1 68 
cells prior to selection at 37°C on LB rich medium con- 
taining chloramphenicol (5 |.ig/ml). Chromosomal DNA was 
extracted from chloramphenicol-resistant clones and ye/A 



gene deletion was verified by PCR using the upstream 
out5 (5'-GAGTTTGACCTGATCATTGCGGC-3') and 
downstream out3 (5'-GAATCTTTTCAGAATATATTG 
CTGG-3') primers. One of the correct strains was dubbed 
168 AyefA-.-.Cm and was used for further study. The 
B. subtilis strain BFS2838 carrying an inactivated trmFO 
(gid-O) gene was kindly provided by S. Serror (European 
functional analysis project of B. .subtilis; http://genome 
.jouy.inra.fr/cgi-bin/micado/fonct_analysis.cgi). 

Escherichia coli strains used in this work are derivatives 
of BW 25113 [¥-, A(araD-araB)567, AlacZ4787{::rmB-3), 
1~, rph-1, A(rhaD-rhaB)568, hsdR514] and were obtained 
from the Coh Genetic Stock Center. The BW 25113 
AtrmAj ArlmCj ArlmD strain was obtained by sequential 
Pl-mediated transductions of the AtrmA753::kan, 
ArlmD783::kan and ArlmC760::kan alleles (16). After 
each transduction, the kanamycin cassette was excised by 
transient expression of the specific FLP recombinase (17). 
Deletion of the three genes was verified by PCR using 
primers adjacent to the inactivated sequences. 

Plasmids 

The 1.4-kb yefA ORE was amplified from genomic DNA 
with the PCR primers atg_Hind (5'CCCAAGCTTATGA 
AAATGAAACCCCCAGTAG-3') and stop_Xba (5'-CT 
AGTCTAGATTATTCTTTTAATTTTATCAACACAC 
AG-3') containing Hindlll and Xbal sites, respectively. 
The PCR fragment was cloned into the Smal-digested 
pUC19 vector to form pUC-je/^ and was then digested 
with Hindlll and Xbal before cloning the yefA sequence 
into these sites in plasmid pDG148 (18) forming pDG- 
yefA. In this vector, expression of yefA was under control 
of the P.spac promoter, and was induced with 0.5-1 mM 
IPTG. The pDG148 and ^BG-yefA vectors were shuttled 
into competent 168 Aj'e/^::Cni ceUs, and transformants 
were selected on LB plates containing chloramphenicol 
(5|ig/ml) and kanamycin (6^g/ml). 

The pBAD-ve/^ vector was obtained by digesting pUC- 
yefA with EcoRI and Xbal, and cloning the yefA frag- 
ment into the same sites in the pBAD24 vector bringing 
yefA expression under the control of the arabinose 
promoter, which was induced with 0.2% arabinose. 
Electrocompetent AtrmA/ArlmC/ArlmD cells were trans- 
formed with the pBAD-yefA vector and selected on LB 
plates containing ampicillin (100 ^g/ml). All recombinants 
were verified by sequencing. 

Preparation of B. subtilis and E. coli rRNA 

The E. coli and B. .subtilis strains were grown at 37°C to an 
ODfioo of 0.6 in 200 ml LB medium containing the appro- 
priate antibiotics. For strains complemented with recom- 
binant copies of yefA and yfjO, IPTG was added at this 
stage to a final concentration of 1 mM to induce protein 
expression. Cells were kept at 37°C for three more hours, 
and were then harvested by centrifugation, washed and 
resuspended on ice in 3 ml of 50 mM Tris-HCl pH 7.5, 
lOmM MgCl2, lOOmM NH4CI before being lysed by son- 
ication. Cell debris was removed by centrifugation and 
ribosomal particles were pelleted by ultracentrifugation 
(19) leaving bulk tRNA in the supernatant. The tRNA 
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and rRNA fractions were extracted with phenol and 
chloroform, the RNAs were recovered by ethanol precipi- 
tation and dissolved in H20. 



MALDI mass spectrometry analyses of rRNA 

RNA was extracted from the ribosomal particles of the 
B. suhtilis 168 wild-type, yefA and yjjO knockouts, and 
from the knockouts complemented with active copies 
of yefA and yfjO; rRNAs were additionally extracted 
from the E. coli ArlmC/ArlmD/AtrmA triple knockout 
and from this strain expressing an active copy of yefA or 
yfjO. The sequence from G725 to G772 within domain II 
of the B. subtilis 23S rRNA (Figure lA) was isolated 
by hybridization to a complementary 48-mer deoxy- 
ohgonucleotide, 5'-CCCACACCTCATCCCCGCACTTT 
TCAACGTGCGTGGGTTCGGGCCTCC, while the 
sequence from C1914 to CI 961 within domain IV of the 
B. subtilis 23S rRNA (Figure IB) was hybridized to the 
48-mer, 5'GTCGGAACTTACCCGACAAGGAATTTC 
GCTACCTTAGGACCGTTATAG. The corresponding 
regions in E. coli rRNA were isolated using similar oligo- 
nucleotides adjusted for the sequence variations shown 
in Figure 1. In each case, lOOpmol of total rRNA was 
hybridized to SOOpmol of deoxyohgonucleotide; unpro- 
tected regions of the rRNAs were digested away with 
mung bean nuclease and RNase A; and the rRNA se- 
quences paired to the deoxyohgonucleotide were sep- 
arated by gel electrophoresis and extracted (20,21). The 
rRNA sequences of around 48 nt were digested with 
either RNase A or RNase Tl in aqueous solution contain- 
ing 3-hydroxypicolinic acid and analysed by MALDI-MS 
(Voyager Elite, Perseptive Biosystems) recording in re- 
flector and positive ion mode (22). Spectra were 
analysed using the program m/z (Proteometrics Inc). 
Tandem mass spectra were recorded in positive ion 
mode on a MicroMass MALDI Q-TOF Ultima mass 
spectrometer (23). 

Analyses of RNAs by HPLC and ion-trap electrospray 
mass spectrometry 

The supernatant fractions containing tRNAs and other 
soluble RNAs were passed through Microcon YM-100 
columns (Millipore) to reduce the amount of impurities. 
The supernatant samples were digested to completion to 
form nucleosides (24), and the pellet samples containing 
the rRNAs were similarly treated. The resultant nucleo- 
sides were analysed using a HPLC/ES-MS ion-trap 
system. Briefly, the loading pump of an Agilent 1100 
HPLC was run at 5 |il/min; the column was a Hypercarb 
150 X 0.3 mm, 5 ^m PGC with a pore size of 250 (Thermo 
Scientific), hnked to an Agilent XCT Ultra 6340 ion trap. 
Separation was achieved on gradients of 0-90% acetoni- 
trile in 0.1% formic acid. ES-MS spectra were recorded in 
negative ion mode scanning initially in the mjz range 220- 
400 for intact nucleosides and then in lower mjz ranges 
after fragmentation (24). Nucleosides from the tRNAs 
were additionally analysed using reversed phase HPLC 
as previously described (24). 
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Figure 1. The three characterized sites of m U modification in 
B. suhtilis and E. coli RNAs. The 23S rRNA secondary structures 
are based on the B. suhtilis sequence and nucleotide differences in 
E. coli are indicated by the arrows pointing to the itahcized nucleotides. 
(A) Portion of domain II of 23 S rRNA, containing hairpins 34 and 
35 and part of helix 33, showing the location of m^U747 modified 
by YefA in B. suhtilis rRNA and RlmC in E. coli. The m'G745 
inodification is seen in E. coli rRNA but is absent in B. suhtilis; the 
presence of pseudouridine >I'746 could not be tested by mass spectrom- 
etry. (B) Region from domain IV of 23S rRNA containing m^U 
1939 that is modified by YefA in B. suhtilis rRNA, and by RlmD 
in E. coli. The B. suhtilis rRNA is methylated on the ribose of nucleo- 
tide CI 920 presumably by an orthologue of the enzyme TlyA (42) and 
this modification is missing in E. coli; the m^C1962 modification is 
present in the E. coli rRNA (47) but was not detected in B. .suhtilis; 
the pseudouridines *1911, *1915 and *1917 could not be verified 
by mass spectrometry. (C) Schematic consensus of tRNA structures 
indicating the U54 target of the E. coli TrmA (48) and the structurally 
unrelated enzyme TrmFO that modifies the same nucleotide in 
B. suhtilis using a tetrahydrofolate cofactor as the methyl donor 
(13,25). 
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RESULTS 

Nucleotide modifications m^U747 and ni^U1939 are 
present in B. suhtilis 23S rRNA 

Isolation of the B. subtilis 23S rRNA domain II region 
(Figure lA) followed by RNase Tl digestion gives rise 
to a unique UUGp fragment containing nucleotide 
U747. In the MALDI-MS spectra of the rRNA from the 



wild-type B. subtilis strain, this fragment produced peaks 
at mjz 990 and mjz 972, respectively, for the hnear and 
cyclic phosphate forms of UUGp (Figure 2A), corres- 
ponding to the trimer plus 14 Da extra mass of one 
methyl group. There was no peak at mjz 976 (or at mjz 
958 for the cyclic phosphate form), indicating that methy- 
lation of the sequence was stoichiometric. Tandem MS 
of this fragment and a longer partially digested RNase 
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Figure 2. MALDI-MS analyses of m U methylation sites in 23S rRNA. The upper row of spectra were generated from 23S rRNA of tlie B. suhtilis 
wild-type strain 168 and show sequences containing m^U nucleotides (shaded in grey) at (A) U747 (in an RNase Tl oligo), (B) U1939 (RNase Tl 
oUgo) and (C) U1939 (RNase A oligo). The corresponding spectral regions from the B. suhtilis AyefAwCm strain are shown in panels (D-F). The 
rRNA spectra from the B. suhtilis AyfjO strain (data not shown) were identical to the wild-type rRNA. The rRNA analysed in spectra (G and H) 
was from B. suhtilis AyefAwCm strain complemented with an active copy of the yeJA gene. (I and J) are spectra from the same rRNA regions of the 
E. coli AiiniCIAiiinD/ AtrmA triple mutant expressing an active copy o{ yeJA (note the E. coli m'G745 is resistant to RNase Tl digestion, producing 
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(>p) and a linear 3'-phospliate (p), the latter having an 18 Da larger mass; the theoretical monoisotopic mass/charge values (injz) are given for the 
predominant peaks. 
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Tl fragment (CACGUUGp) showed that the methyl 
group resided on nucleotide U747 (Supplementary 
Figure SI). 

In domain IV of the rRNA from the wild-type 
B. subtilis strain (Figure IB), RNase Tl digestion 
produced a unique decamer AAAUUCCUUGp contain- 
ing nucleotide U1939. The MS peak at mjz 3199 (and at 
mjz 3181 with a cyclic phosphate) corresponds to the 
decamer plus a single methyl group (Figure 2B). The 
lack of any peak at mjz 3185 (and at mjz 3167 with the 
cyclic phosphate) indicated stoichiometric methylation of 
this sequence. The RNase A fragment GAAAUp (nucleo- 
tides 1935-1939) flew at mjz 1671 (Figure 2C) showing 
that this overlapping sequence also carried the methyl 
group. The exact location of the methylation site was pin- 
pointed to nucleotide U1939 by tandem MS analysis of 
this sequence using the sUghtly longer fragment GAAAU 
Up formed by partial RNase A digestion (Supplementary 
Figure S2). 

These regions of the B. subtilis rRNA were analysed by 
primer extension, which did not reveal any reverse tran- 
scriptase stops or pauses at nucleotides U747 and U1939 
(data not shown), and ruled out that the metliylations 
could be at the N3 of the bases or the 2'-OH of the 
riboses (21). These findings were fully consistent with 



the location of methyl groups added by YefA being at 
the C5-positions of nucleotides U747 and U1939. 

Candidate methyltransferases for the m^U747 and 
m^U1939 modifications 

A BLAST search of the B. subtilis genome using the 
sequence of the tRNA-modifying enzyme TrmFO as 
the query revealed no paralogue that might function in 
rRNA modification (25). Similar searches using the 
E. coli 23S rRNA m^U methyltransferases RlmC and 
RlmD, revealed two B. subtilis orthologues, YefA and 
YfjO. These two putative B. subtilis methyltransferases are 
43% identical (68% similar) in their amino acid sequences 
(Supplementary Figure S3) and are clearly paralogues. 
Comparison with the E. coli methyltransferases revealed 
that YefA displays 22% and 30% identity (40% and 49% 
similarity) to RlmC and RlmD, respectively (Figure 3); 
YjfO is slightly less close to the E. coli se- 
quences and shows corresponding identities of 19% and 
24% (42% and 47% similarity). Phylogenetically, aU the 
four proteins YefA, YfjO, RlmC and RlmD fit into the 
same cluster of orthologous groups COG2265 (11). The 
yefA and yfjO genes were therefore candidates (and, in 
fact, the only likely candidates) for encoding enzymes 
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Figure 3. Amino acid sequence alignment of m U methyltransferases. The E. coli RlmC and RlmD enzymes, which are respectively specific for 23S 
rRNA nucleotides U747 and U1939, are aligned with the B. subtilis YefA sequence. Identical amino acid residues are highlighted in black, and grey 
indicates highly similar residues. Full-length sequences are shown; RlmC lacks the N-terminal extension seen in RlmD and YefA. The common core 
sequences were used for calculating percentages of amino acid identity and similarity. 



Nucleic Acids Research, 2011, Vol. 39, No. 21 9373 



that could catalyse the m^U modifications in B. subtilis 
23S rRNA. 

Identifying the enzyme(s) responsible for the B. subtilis 
ni^U747 and ni^U1939 modifications 

Strains BG12826 (ye/y4::pMutin2) and BG12911 {yfjO:: 
pMutin2) were procured and then tested by PGR to estab- 
lished whether gene inactivation had been successful. The 
BG12911 0f//'C::pMutin2) strain produced PGR fragment 
patterns that were fully consistent with targeted inactiva- 
tion of the yfjO gene; however, the (j'<?/^::pMutin2) strain 
showed anomalies indicating that genes neighbouring 
yefA on the chromosome might also have been disturbed. 
We therefore made another yefA knockout (Aj'e/^::Gm) 
in-house using a different procedure; subsequent PGR 
analysis (data not shown) confirmed that yefA had been 
inactivated without disruption of any other gene. 

Analysis of the rRNA from the l^yefAwCm strain by 
MALDI-MS showed that methylation had been lost at 
both U747 (Figure 2D) and U1939 (Figure 2E and F). 
Parenthetically, we note that the initial yefA knockout, 
BG12826 (>'e/L4::pMutin2), had also lost the methylations 
at both of these nucleotides (data not shown); however, we 
continued all recombinant work with our own Aj'e/^::Gm 
strain. Transforming this latter strain with an active copy 
of yefA on a plasmid restored to a significant degree 
the' methylation at both U747 (Figure 2G) and U1939 
(Figure 2H). 

The dual methylation function of YefA was addition- 
ally tested in a strain of E. coli in which all three m^U 
methyltransferase genes had been deleted. In this heterol- 
ogous expression system, methylation was partially 
restored at both U747 (Figure 21) and U1939 positions 
in the E. coli 23S rRNA (Figure 2J) with the latter site 
being more efficiently methylated. Methylated uridines 
from E. coli and B. subtilis tRNAs were analysed using 
an HPLG/ES-MS ion-trap system, which verified that the 
methyl group was on the G5-position of the uracil base 
(Supplementary Figure S4). Further HPLG analysis 
showed that YefA is not involved in m^U modification 
of tRNAs (Supplementary Figure S5). 

Parallel sets of experiments were carried out in attempts 
to identify the function of the YfjO enzyme. In the 
BG12911 strain, in which this enzyme is missing, nucleo- 
tides m^U747 and m^U1939 remained fully methylated. 
Furthermore, HPLG analyses of tRNAs from the 
B. subtilis strains showed that, similar to YefA, YfjO plays 
no role in tRNA m^U modification (Supplementary Figure 
S5). At present therefore, the methylation target for the 
YfjO enzyme remains unknown. 

DISCUSSION 

The rRNAs of all Hving organisms are post- 
transcriptionally modified to improve their performance 
in protein synthesis (2,26), although there are marked dif- 
ferences in the mechanisms by which rRNA modification 
is achieved. In Archaea and Eukaryota, most rRNA modi- 
fications are uridine isonierizations and 2'-0-ribose 
methylations (8). These modifications are generally 



added by enzyme complexes that are guided to their nu- 
cleotide targets by complementary sequences in the 
sRNAs of Archaea (27,28) and the snoRNAs of 
Eukaryota (29-31), although some non-guided modifica- 
tions have also been noted (32,33). In the rRNAs of 
Bacteria, uridine isomerizations and 2'-(?-methylations 
also occur, although base methylation is the most 
frequent type of modification (4,34,35). None of the bac- 
terial rRNA modifications requires the help of guide 
RNAs, and the general view has been that each modifica- 
tion is added by a single-specific enzyme that is capable of 
independently finding its target nucleotide. 

On the whole, this view has been substantiated by 
studies of the numerous rRNA modifying enzymes in 
E. coli [reviewed in (34)] with only a few amendments 
being necessary. For instance, a handful of enzymes inclu- 
ding the pseudouridine synthase RluD (20,36), the highly 
conserved m^A dimethyltransferase RsmA/KsgA (37-39) 
and the m^G methyltransferase RsmF of Thermits 
thermophilus (40) modify multiple sites that are immedi- 
ately adjacent (RsmA) or within several nucleotides of 
each other on the bacterial rRNAs (RluD and RsmF). It 
could be argued that these enzymes engage in a single 
binding event with the rRNA during which a second or 
third nucleotide is accommodated into the active site 
without the need for the enzyme to dissociate from the 
substrate. Rare exceptions include the 2'-0- 
methyltransferase TlyA from Mycobacterium tuberculosis, 
which in separate reactions modifies nucleotides on both 
subunits of the ribosome (41), RluG which catalyses 
pseudouridine formation at distant locations in E. coli 
23S rRNA (42), and RluA which forms pseudouridines 
in both the tRNAs and 23S rRNA of E. coli (43). 

The data presented here show that B. subtilis YefA is 
able to recognize, bind and modify two targets that are at 
separate locations. From the MS analyses, it can be seen 
that B. subtilis 23 S rRNA possesses the modifications m^U 
747 (Figure 2A) and m^U1939 (Figure 2B and G) and that 
methylation is stoichiometric at each of the two 
nucleotides. Both methylations are lost upon inactivation 
of yefA and both are then restored upon complementation 
with an active copy of this gene. However, as it could 
formally have been argued that recognition by YefA of 
the two sites was facilitated by an auxiliary B. subtilis 
factor or was perhaps due to some peculiarity in the B. 
subtilis rRNA, a strain of E. coli was constructed in which 
all three of its m^U methyltransferase genes (rlmC. rlmD 
and trmA) were deleted. Upon heterologous expression of 
yefA in this E. coli strain, methylation at U747 and U1939 
was regained showing that the specific recognition of 
dual targets by YefA is a property inherent within the 
structure of this enzyme. Analyses of the tRNA compo- 
nent of these cells showed that YefA is not involved in 
tRNA methylation. Thus, in B. subtilis, YefA is specific 
for the 23S rRNA modifications m^U747 and m^U1939, 
while modification at m^U54 tRNAs is added exclusively 
by the folate-dependent methyltransferase TrmFO (13,25). 

In wild-type E. coli, the G5 positions of nucleotides 
U747 and U1939 are stoichiometrically methylated by 
the independent action of the enzymes RlmG and RlmD 
(9,10). The requirement for two autonomous m^U rRNA 
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methyltransferases in E. coli, while a single enzyme can do 
the same job in B. suhtilis, raises the question as to how 
the functions of YefA, RlmC and RlmD evolved. One 
likely scenario, consistent with phylogenetic analyses 
(3,5), would be that the contemporary B. suhtilis YefA 
enzyme arose from a more promiscuous multi-site, 
AdoMet-dependent m^U rRNA methyltransferase that 
was an earlier (and now perhaps extinct) member of the 
COG2265 cluster of orthologous genes. Similarly, in E. 
coli and related species, single-site specificity evolved fol- 
lowing duplication of an ancestral COG2265 gene giving 
rise to the rlmC and limD paralogues. The limited amino 
acid identity between the three enzyme sequences (Figure 
3) favours the idea that separation of RlmC, RlmD and 
YefA was indeed an ancient event. In the context of func- 
tional shifts, we note that while YefA stoichiometrically 
methylates both of its B. sublilis XilAl and U1939 targets 
(Figure 2A and B) and the recombinant YefA methylates 
with high efficiency at U1939 in B. suhtilis (Figure 2H) and 
E. coli (Figure 2J), the enzyme was clearly less effective at 
recognizing U747 in E. coli where an appreciable residue 
of unmethylated nucleotide remained (Figure 21). A pre- 
cedence for target shifts in m^U RNA methyltransferases 
has already been observed in the archaeon Pyrococciis 
ahyssi where two functional RlmD-like paralogues are 
present, and neither of these modified U1939 (3,5). 

Although m^U RNA methyltransferases have changed 
their modification targets during evolution, the require- 
ment for m^U modification at specific locations has 
remained remarkably constant. The majority of organisms 
appear to possess a homologue of either TrmA or TrmFO 
for m^U54 modification in tRNA (25). Likewise, the 
majority of bacteria encoded at least one RlmD gene, 
which fits with the observation that all bacteria studied 
to date possess the 23S rRNA m^U1939 methylation 
and many also have the m^U747 modification. The con- 
servation of modified nucleotides at identical positions in 
rRNA would suggest that they play essential roles. 
Nevertheless, the E. coli knockout strain grew comfort- 
ably (in rich medium) despite the loss of all three of its 
m^U's and, even though the strain would be expected to 
be at a disadvantage if subjected to stresses including anti- 
biotics (44), its resilience was noteworthy. Similarly, the 
deletion of yefA (or yfjO) in B. suhtilis had no major effect 
on the viability of the strains under standard laboratory 
growth conditions. 

A functional role for the paralogue YfjO remains to be 
determined. The YfjO enzyme does not play a role in m^U 
methylation of B. suhtilis tRNAs (Supplementary Figure 
S4). The HPLC and ion-trap-MS approaches 
(Supplementary Figures S4 and S5) are less well suited 
to the analysis of m^U nucleotides in the larger rRNAs, 
where this modification makes up a much smaller a pro- 
portion of the total uridine complement and, although we 
did show by MALDI-MS that YfjO does not modify 
within the 747 or 1939 regions of 23 S rRNA, we have 
not ruled out that this enzyme might have an rRNA mod- 
ifying function. In fact, four rlmD orthologues and one 
trmFO orthologue are evident in the genome of 
Aclioleplasma laidlawii (45), which belongs to the 
Mollicutes and is distantly related to B. suhtilis, and this 



fits well with an earlier analysis that indicated up to six 
m^U nucleotides at unidentified locations in the A. 
laidlawii 23S rRNA (46). The present evidence is consist- 
ent with certain members of the Gram-positive Firniicutes 
(which include Bacillus and Acholeplasma) having, in 
addition to U747 and U1939, m^U modification(s) in at 
least one other region of the 23 S rRNA. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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